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We report on the status of the determination of the heavy quark masses from our calculation of the quarkonia 
spectra. All sources of systematic errors that enter the quark mass determination are accounted for. We explicitly 
keep ma 7^ in the perturbative calculation relating the bare lattice mass to a renormalized mass. Our results 
are still preliminary. 



1. INTRODUCTION 

Quarkonia, mesons containing a heavy quark 
and anti-quark, are at present the best under- 
stood hadronic systems. As has been argued by 
Lepage [1], quarkonia are also the easiest systems 
to study with lattice QCD, and systematic errors 
can be analyzed using potential models. Con- 
trol over systematic errors in turn allows the ex- 
traction of Standard Model parameters from the 
quarkonia spectra. 

We report on the status of the charm-quark 
mass determination from the cc spectrum using 
the Fermilab action [2]. This action is a gener- 
alization of previous approaches, which encom- 
passes the non-relativistic limit for heavy quarks 
as well as Wilson's relativistic action for light 
quarks. Lattice-spacing artifacts are analyzed for 
quarks with arbitrary mass. The determination 
of the strong coupling from the cc and hh spectra 
has already been presented elsewhere [3]. 

2. THE SPECTRUM 

Figure 1 shows the cc spectrum calculated with 
the Fermilab action in the quenched approxima- 
tion [4]. The lattice spacing is set with the hg—lS 
splitting (see Refs. [3,4] for details on the numer- 
ical calculations). 

Let E denote the energy of the quarkonium 
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Figure 1. The charmonium spectrum in compar- 
ison with experiment. 

state, extracted from the exponential fall-off. The 
rest mass. Mi = E{p = 0), and the kinetic mass, 
g, are unequal [2] for the 
actions used in our numerical calculations. This is 
acceptable in non-relativistic systems if the bare 
lattice mass is tuned such that the kinetic mass, 
M2, equals the desired ground-state mass. Fig- 



Table 1 

Comparison of higher order lattice spacing errors 
for the IS state and the IP-IS splitting. 



(GeV) 


c 


IS (MeV) 


IP-IS (MeV) 


2.4 


1.4 


-F30 


- 4 




0.0 


-F50 


-30 


1.8 


1.4 


-F40 


- 7 




0.0 


-F70 


-40 



Q I I I I I I I I I I I 

1 2 

Moa 

Figure 2. Numerical results for Mia (+) and 
(□) vs. Moa. 

ure 2 shows the masses Mi and M2 as functions 

of the bare lattice mass Mq = 7^ in lattice 

spacing units for the 24^ X 48 lattice at /? = 6 . 1 [4] . 
The effect of higher order lattice spacing errors 
on the spectrum can be estimated using poten- 
tial models. Table 1 compares the effect for the 
spin-averaged IS state and the IP-IS splitting in 
the cc spectrum at several lattice spacings and for 
the Wilson and the improved (c = 1.4) action [4]. 
Richardson potential model wave functions were 
used for the estimate. 

3. THE SELF ENERGY 

The calculation of the self energy with nioa ^ 
has been outlined in [5], where preliminay results 
at one-loop for the Wilson action were shown. 
In order to distinguish quark and meson masses, 
we shall denote meson masses by Mi 2 and quark 
masses by mi 2- The perturbative expansion for 
mi is written as 

mia = \og{\^Moa)^g'^m^^\^0{g'^) . (1) 

Figure 3 shows the one-loop term, rn:^\ for the 
Wilson and improved (c = 1.4) actions as a func- 
tion of rn'f ^ = log(l-l-Moa) with and with- 
out tadpole improvement (using the plaquette) 
[8]. The corresponding effective scales for the 
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Figure 3. 'm'j^\ as function of m'^'^-'a for the Wil- 
son action (solid), improved action with c = 1.4 
(long dashed) and with tadpole improvement for 
both actions (dash-dotted, short dashed). 

coupling, q* , as defined in Ref. [6] (eq. (19)) are 
shown in Figure 4. 

4. RESULTS 

The determination of the charm quark mass 
proceeds analogously to the h quark mass deter- 
mination of Ref. [7]. The binding energy is ob- 
tained from -Ebind = Ml — 2mi , which together 
with the experimentally observed spin-averaged 
mass of the IS state, Mj'g'' = 1/4(M^^ 3Mj/^), 
determines the charm quark pole mass as 

<°'^ = ^(Ks'^-^H„d) . (2) 

Alternatively, after tuning the bare quark mass, 
such that M2 = Mj^g'^, the pole mass can be ob- 
tained from: 

"^ch'*" = "^2 = Zmm'^°^ , (3) 

with Zfn = 1 + Zrn^ at one-loop. Preliminary 
results for Zm^ exist. Once the calculation is com- 
pleted [8], the comparison of eqs. (2) and (3) will 
give an important consistency check. 

Of course, as always, all systematic errors aris- 
ing from the lattice QCD calculation need to be 
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Figure 4. q* for m^^^ as function of m^^^ a for the 
Wilson action (solid), improved action with c = 
1.4 (long dashed) and with tadpole improvement 
for both actions (dash-dotted, short dashed). 

under control for a phenomenologically interest- 
ing result. In particular, the systematic error in- 
troduced by the omission of sea quarks has to 
be removed. The short-distance corrections that 
introduced the dominant uncertainty to the 
determination from quarkonia [3] are absent for 
the pole mass determination, because this effec- 
tive mass does not run for momenta below it's 
mass. This argument has been tested from first 
principles in Ref. [7] for the h quark pole mass 
with consistent results. 

The statistical errors are dominated by the 
IP-IS splitting which determines a~^, and are 
< 10 % of the binding energy. An estimate of the 
higher order lattice spacing errors to the spin- 
averaged spectrum is listed in table 1. Finite vol- 
ume errors are negligible. 

The dominant systematic error comes from the 
unknown higher order corrections to eq. (1) (and 
eq. (3)); a preliminary estimate is 100 — 200 MeV. 
Our preliminary result is m^^'^ = 1.5(2) GeV. 

The MS mass for the charm quark has been de- 
termined from a compilation of D meson calcula- 
tions in the quenched approximation [9] . However 
not all systematic errors were considered. In par- 
ticular, sea-quark effects cannot, in this case, be 
estimated phenomenologically, leaving this sys- 



tematic error uncontrolled. Once our analysis for 
the MS mass is completed, we will be able to com- 
pare our result with that of Ref. [9]. 

The analysis of the hh spectrum for a determi- 
nation of the h quark pole (and MS) mass is also 
in progress. 

As already pointed out in Ref. [7], the heavy 
quark mass determination is limited by our lack of 
knowledge of the higher order perturbative terms. 
The extension of the self energy calculation to 
two-loops will require tools different from those 
used in the one-loop calculation [10]. 
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